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of the apical oxygen from O to O should lead to a reduced phonon frequency 12, and hence 
to an increased renormalization factor in which D. enters exponentially. In this case a shift of 
the main peak to higher frequency would be expected. Furthermore, one should observe a sign 
of short-range CDW order resulting from the polaron-polaron interactions in this case. 
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open boundary conditions, which is the simplest way to include the effect of the oxygen va- 
cancies. We found that the effective electron correlations lead to the appearance of the high- 
frequency tail also in this situation, provided the density is close to 1/2, which represents the 
estimated experimental value. In the low-frequency region, we also found the characteristic 
peak, the origin of which can be traced back to the lowest particle-hole excitation of the non- 
interacting system. At even lower energies, the absorption due to the renormaUzed motion of 
the domain walls partly closes the quasi-gap. 

When faced with the experimental spectrum, the agreement is fairly good. While the 
high-frequency tail is reproduced very well, at low-energy there is some discrepancy: The 
conductivity falls off two sharply on the low-energy side of the main peak, and the quasi- 
gap leaves some missing weight. Most likely, this shortcoming is due to the fact that only 
chains with length N<14 could be considered in our calculations. The contribution of longer 
chains has to be expected to close the quasi-gap, since its size diminishes for increasing chain 
lengths. Adding the spin degrees of freedom was shown to lead to additional absorption in the 
low-frequency region, and a slight shift of the main peak to lower energies, without strongly 
affecting the high-frequency tail. However, for physical values of J/t < 0.5, the overall in- 
fluence of the spin degrees of freedom was found not to be very significant, as expected due 
to charge-spin separation. The inclusion of additional phonon states mainly has the effect of 
renormaUzing the e-ph coupUng strength to larger values. Apart from a sUght smoothening of 
the spectra, it does not alter the results. 

In conclusion, there are two likely mechanisms which could give rise to the observed 
conductivity of the CuOj -chains in the 1-2-3 materials: While the first one (strong potential 
disorder) gives excellent agreement with the experimental spectrum it suffers from the difficult 
justification of the required disorder strength. The proposed polaronic e-ph coupUng can be 
well justified assuming reasonable physical parameters, and also gives fair agreement with the 
experiment. 

It would be highly desirable to have an experimental verification for one or the other mech- 
anism. The total spectral weight expressed in the effective carrier number per unit cell was 
found to be very similar in both models n^ff = 0.4-0.45, and compared well with experiment 
n^ff -0.5. One possibiUty to distinguish the two mechanisms could be the investigation of the 
'isotope effect' on the conductivity spectrum. If the e-ph model is appUcable, the substitution 
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show (i) a low-energy peak, the position of which depends on the disorder strength, and (ii) a 
high-energy tail. Therefore it appears to be a promising candidate to explain the experiment, 
hideed, we found that it is possible to fit the measured spectrum very accurately by this model. 
However, a quantitative analysis revealed that a very large disorder strength (W ~ 1.2- 1.6eV) 
is needed to obtain the correct peak position. There is no experimental evidence for such strong 
disorder in the real materials, so it remains unclear whether this is the correct model. A possible 
type of defect, which could in principle produce such strong potentials, is an interstitial oxygen 
atom between two Cu(l) belonging to neighbouring chains. But the density of such defects 
(= 10%) required to produce such a strong random potential on each site seems unphysically 
large. 

As a third likely source of strong scattering in the chains we considered polaronic electron- 
phonon coupUng. In particular, we argued that the phonon mode which involves the displace- 
ment of the apical oxygen (0(4) site) in c-direction, should couple rather strongly to the holes 
in the chains. We showed that this coupUng is well described by an additional Holstein term 
in the Hamiltonian, the local ion displacement being linearly coupled to the local charge den- 
sity. For simplicity, we assumed this phonon mode to be dispersion-less. Again neglecting the 
spin degrees of freedom, the effective model becomes the standard Holstein molecular-crystal 
model for spinless fermions. 

In consideration of the importance of the effective electron-electron interactions, we chose 
a numerical technique for the calculation of (7(<o), which is capable of deaUng with such 
interactions in a clean way, the Lanczos-method. Apart from the limitation to finite lattices, the 
drawback of this method when applied to an e-ph Hamiltonian is, that one is forced to restrict 
the possible phonon states to some finite subspace, in order to obtain a finite-dimensional total 
Hilbert space. 

In an accompanying paper [25], we investigated different possible variational subspaces 
for the phonons, each of them allowing two states at a given lattice site. The best choice for 
the two local phonon states, which we called Ansatz II, was found to be the vacuum states of 

the undisplaced, and the displaced oscillator |0)/, |0)/. The agreement from Ansatz 11 of the 
parameters entering the effective Hamiltonian with the exact ones, as well as with the QMC 
data was excellent. 

For a comparison with the experimental spectrum, we used Ansatz II in the presence of 
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VII. Conclusions 



In this article, we discussed different ptiysical mechanisms which could be responsible for the 
experimentally observed frequency dependence of the optical conductivity in the CUO3 -chains 
of the 1-2-3 compounds. The main features of the conductivity are a broad mid-infrared peak 
at (0 = 0.2eV, and a subsequent slowly falUng (much slower than l/o ) high-frequency tail. 

In view of the similarity between the CuO bonding patterns of the CuOj -planes, and the 
CuOg-chains, we argued that a ID ?-/-model, derived from a multi-band Hubbard model should 
capture the relevant low-energy physics of the chains. However, the known fact that this model 
does not show any substantial optical absorption (due to the charge-spin decoupling in ID), 
forced us to consider additional terms in the Hamiltonian, which could lead to the observed 
strong scattering, and the resulting absorption. 

The inherent tendency of oxides towards the formation of defects naturally suggests to 
include the resulting perturbations of the charge motion in a model of the electronic structure. 
In the case of the 1-2-3 compounds, the oxygen vacancies on the 0(l)-site connecting two 
neighbouring Cu(l) in the chains present the most serious and also most common type of 
defect. 

The simplest way to model the vacancy disorder is to assume that a vacancy acts as an 
almost impenetrable barrier. This leads to a decomposition of an individual chain into effec- 
tively discoimected segments, between which charge motion is severely hindered. This effect 
can be incorporated in the Hamiltonian by a strongly reduced hopping matrix element between 
the lattice sites next to a vacancy. For the optical conductivity, this has the consequence that 
the dc conductivity vanishes, and its weight is transferred to finite frequencies. 

Neglecting the magnetic term in the ?-/-model, we calculated the optical conductivity for 
such a model, assuming a uniform random distribution of vacancies, characterized by the mean 
length of an individual segment. Fixing the energy scale by setting t = 0.4eV, we found that 
this model exhibits a low-energy peak, but it is neither able to correctly reproduce the peak 
position nor does it show the observed high-frequency tail. For a more realistic model, where 
in addition to the disorder in t, we also included a repulsive potential on the sites neighbouring 
the vacancy, we did not find any better agreement with the experimental spectrum. 

As a second possibiUty, we included potential disorder, which can be modeled by random 
on-site energies on each lattice site. The ac conductivity of such models is well known to 
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Figure 18: A comparison of Oi (©) obtained from Ansatz n at different values of the spin exchange 
coupling J = 0.0, 0.5, 1.0, and fixed N = 12, n,, = 5, Eg = 0.65, Q. = 0.2, = 111 using OBC. The inset 
contrasts the case of 7 = 0.5 with a calculation for y = 0, £5 =0.55. 



37 



showed in Ref.[25], at small density, and QJt < 0. 1, excitations on the scale Q. become relevant, 
and this leads to a qualitative difference between the results from Ansatz II and those from the 
true Hamiltonian. Note, also, that the additional phonon states lead to a slight smoothening of 
the spectra as compared to Ansatz II. 

(ill) Finally, we should consider the possibility that the scattering from spin exchange could 
be enhanced by the e-ph coupUng as compared to the simple t-J model (see the discussion 
in section II). Naively, one expects a competition between the spin and the polaronic e-ph 
coupUng. The former favours singlet bonds between electrons on neighbouring sites, i. e., an 
effective attraction, whereas the latter leads to an effective nearest neighbour repulsion, i. e., 
CDW correlations. This competition might lead to interesting effects, which we investigate in 
more detail in Ref. [25]. 

In Fig. 18 we show the spectra calculated for a 12-site lattice with - 5, total spin 
Z-component of - 1/2, parameters Eg - 0.65, i2 = 0.2, and different values of the spin 
exchange coupUng /. The ground state for all the parameters has the lowest possible value 
for the total spin, S = 1/2, as found in the simple t-J model. We note that the spin degrees of 
freedom enhance the low-frequency weight below the main peak, and shift the latter to lower 
energy, while the high-frequency weight is reduced. The additional low-frequency weight in- 
creases as a function of /. This effect of the spin exchange on the low-frequency region is 
very similar to what happens when the e-ph coupling strength is reduced. A possible explana- 
tion for this behaviour could therefore be that the nearest neighbour spin exchange coupling, 
which favours electrons on neighbouring sites, leads to an effective weakening of the repul- 
sive CDW-correlations on the scale <o < /. This can be verified in the inset of Fig. 18, where 
one observes a strong similarity between a spectrum at / # and some E^, with a spectrum at 
J = and a smaller E'^ < Eg. Apart from this effect, the spin degrees of freedom do not play 
a significant role. 

From these considerations, it seems most likely that the missing weight should come from 
the contributions of larger chains. We expect that this should fill in the quasi-gap. Even 
though we cannot expUcitly calculate a spectrum for this situation, our results strongly suggest 
that a spectra obtained from a superposition of random length chain segments (similar to the 
case we calculated in section III without e-ph coupUng) would give good agreement with the 
experimental one. 
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Figure 16: A comparison of spectra al''^{(o) obtained 
for different combinations of lattice sizes and fillings 
[N, He] = ([10, 4], [12, 5], [14, 6]) using OBC and fixed 
values of Eb = 0.65, Q. = 0.2. 
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0.2, and different coupling strengths. The e-ph coupling was chosen such that the qualitative 
shape of the peak and the high-frequency tail are similar for the two approximations. We 
observe an excellent agreement between the two cases over the entire frequency range, i. e., 
the low-frequency weight is not affected by the reduction of the phonon Hilbert space. This 
agreement supports the hypothesis [25] that the neglected phonon states in Ansatz II can be 
taken into account by a renormaUzed e-ph coupUng. In other words, the spectrum we obtain 
from a calculation using Ansatz 11 at coupUng presents an excellent approximation of the 
true spectrum (using the full Hilbert space) at a coupUng strength E'^>Eg. At least this should 
be the case if both, the particle density n, and the phonon frequency Q. are not too small. As we 
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Figure 17: A comparison of Oi^ico) obtained from Ansatz II with the results using a four-phonon basis 
for different couplings. The spectra were obtained for AT = 9 with OBC at = 4 and €1 = 0.2. 
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Figure 15: A comparison of the experimental spec- 
trum obtained in Ref. [1] with two curves calculated 
using OBC and parameters (a) A'^ = 14, Eb = 0.65, and 
(b) N = 13, Eb = 0.55, both at n« = 6, ^ = 0.2. 
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Qq the lattice constant along the chain, i. e., the b-direction. From our calculations Fig. 15, 
we obtain a value of n^^ ~ 0.43 - 0.45. This compares rather well with the experimental one 
n^fi -0.5. The discrepancy comes from the low-frequency part. 

There are three possible explanations for the shortcoming at low frequencies: (i) Since 
in our model, the position of the main peak, and therefore the size of the quasi-gap is deter- 
mined by the energy of the lowest particle-hole excitation of the corresponding non-interacting 
system, it is clear, that it should diminish as a function of increasing system size. Hence, it 
is plausible to assume, that the contribution of longer chains (A^ > 14) to the total spectrum 
in the real samples would fill in the missing spectral weight at low energy. Such long chain 
segments should certainly be present in these materials, provided that the oxygen vacancies 
occupy random positions, and do not form an ordered superlattice, for which there is no ex- 
perimental evidence. Unfortunately, the limitations posed by the available computer power, 
make it impossible to demonstrate this explicitly by a calculation. However, the comparison 
in Fig. 16 between the spectra for increasing system size clearly shows the shift of the peak. 

(ii) The second explanation could be that the approximation of Ansatz 11 sUghtly underes- 
timates the contribution of renormaUzed (co < (0^^^) versus unrenormaUzed (weight of main 
peak) charge motion, so that in reality (using the full phonon Hilbert space) there would not be 
a quasi-gap, but rather a smooth decay of the conductivity as we approach co -0. If this was 
the case, we should see at least the onset of such a shift of spectral weight, when comparing 
the results from Ansatz II with the calculations performed using the two additional excited 
phonon states |1);, |1); . 

Fig. 17 shows a comparison between the two approaches for a 9-site lattice with = 4, H = 
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Figure 14: The frequency dependence of Oi (©) obtained for a lattice with N = 14 using OBC. The 
plots show spectra calculated for different fillings, and coupling strengths Eg (close to at fixed phonon 
frequency i2 = 0.15. 
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where Mq the free electron mass, V^eu the volume of the unit cell, e the elementary charge, and 
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Figure 13: The low-frequency part of CT] (co) ob- 
tained for a lattice with = 12 using OBC at = 5. 
The plots show spectra calculated for different cou- 
pling strengths Eb at fixed phonon frequency H = 0.2. 
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phonon frequency changes the renormahzation factor = e~^^'^. 

The inspection of the spectra at different doping levels showed that the high-frequency 
tail, which is a characteristic feature of the experimental spectrum, is present only for densities 
around 1/2. This is a promising feature of the e-ph model, since we know from experimental 
estimates, that the hole density in the Cu03-chains of both, PrBa2Cu307_g and YBa2Cu307_5 
is approximately 1/2 per Cu, which translates to the value n = 1/2 in our model. Hence, we 
shall now compare the experimental spectrum of PrBa2Cu30y_^ with those calculated spectra 
which give the best agreement. 

Fig. 15 shows the experimental gI^^{(0) together with the ones calculated for parameters 
(a) N =U,n^ = 6,£l = 0.2, = 0.65, and (b) N = 13, = 6, £1 = 0.2, = 0.55. We 
have used the usual value t = 0.4eV to fix the energy scale. The agreement of the calculated 
with the experimental curve is fairly good. In particular, the high-frequency tail is very well 
approximated by the calculation, and also the ratio g(co - (Op^.^]J/a{co = l.OeV) of the peak 
height to the weight at the bottom of the tail (co - l.OeV) is correct. However, there is some 
discrepancy at small frequencies: The quasi-gap below the main peak is not observed in the 
experiment. There, the conductivity drops very sharply when approaching zero frequency. 

The total spectral weight can be expressed in terms of an effective carrier number per unit 
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pearance of the low-energy peak can clearly be related to an optical excitation, which is al- 
ready present in the non-interacting case, namely the lowest particle-hole excitation of the free 
electron system. The corresponding energy is largest at half-filling (due to the simple cosine 
band- structure), and decreases monotonically with the density. This peak dominates the opti- 
cal conductivity in the presence or absence of the e-ph interactions, at least for not too strong 
coupUng (where it is suppressed for densities close to 1/2). However, in the interacting case, 
the peak is shifted to higher energies (as compared to the free system), and also broadened as 
the coupUng increases. 

(ii) The second new feature is the appearance of spectral weight at very low energy below 
the dominant peak, which partly fills in the quasi-gap. This weight decreases as a function of 
increasing coupling in much the same way as the weight under the peak decreases. The com- 
bined reduction of weight of these two features can be related to the corresponding decrease of 
the Drude weight obtained for PBC. Physically, the low-energy weight can be associated with 
the renormaUzed charge motion, which occurs on the energy scale ~ t. For PBC, this motion 
contributes to the dc conductivity, whereas here, due to the reflecting walls, the corresponding 
spectral weight appears at finite frequencies of the order of co ~ t. 

This effect can be observed, when we concentrate on the low-energy spectrum below the 
main peak, and use a smaller broadening 5 = 0.01 for our plots. Fig. 13 shows this part of the 
spectrum for the case N = \2, = 5, Q. = 0.2, and three different coupUng strengths. Note, 
how the dominant peak in this low-energy part is shifted towards lower frequencies as the 
coupling is increased, according to the dependence oft= te~^^'^. A detailed analysis of the 
peak position revealed exactly the exponential dependence. Furthermore, for strong coupling 
Eg > E^g, the inspection of the change in this low-energy spectral weight as a function of doping 
reveals a very similar dependence as obtained for the Drude weight in the case of PBC (see 
Fig. 11): For densities close to 1/2, this weight increases with the number of domain walls, 
whereas for n > 0, it increases with the number of electrons. In summary, we can say that the 
low-energy spectrum clearly shows the presence of simultaneous renormaUzed (<o ~ t), and 
unrenormalized charge motion (co ~ t). 

In Fig. 14, we show the same spectra for lower phonon frequency Q.-0. 15, and various 
coupUng strengths. The quaUtative shape of g{(0) is very similar to the previous case, but note, 
that the coupUng strength needed to obtain the high-frequency tail is now smaller, because the 
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Figure 12: The frequency dependence of Oi (©) obtained for a lattice with N = 14 using OBC. The 
plots show spectra calculated for different fillings, and coupling strengths Eg (close to at fixed phonon 
frequency Q, = 0.2. 



tail follows smoothly onto the main peak, whereas for = 4, 5, the peak looks more like an 
isolated structure superposed with the slowly falling background absorption. 

The differences compared to the spectra with PBC in ref. [25] are two-fold: (i) The ap- 
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VI. Comparison with experiment 



In this section, we restiict ourselves to a set of parameters for the Holstein Hamiltonian, which 
represents those expected in the real chains. In section III, we already pointed out the im- 
portance of defects in a realistic model for the chains. We also mentioned the experimental 
finding that the oxygen vacancies on the 0(l)-site in the chains, are by far the most important 
type of defects. Such a vacancy presents a serious perturbation to the electron motion, since 
it effectively cuts the chains into more or less disconnected segments. In consideration of this 
fact, we think it is justified to assume that the dominant part of the optical response originates 
from the excitations in finite chain segments, and that the inter-segment hopping only plays 
the role of a small perturbation, which we neglect. For our calculations, this means that we 
should use open boundary conditions (OBC), i. e., an oxygen vacancy is modelled by a purely 
reflecting wall. 

We already studied a similar case in section III, where we considered electrons without 
coupling to phonons. One obvious property of electronic models with OBC is that they cannot 
support any dc current. Hence, as already mentioned earlier, the Drude weight is transferred to 
finite frequencies, and especially for rather short chain lengths, this leads to a dramatic change 
in the low-energy response. 

In this section, we shall study the optical conductivity of the Holstein Hamiltonian in the 
presence of OBC. We shall restrict the parameters to be f = 0.4eV, and DJt = 0.15 - 0.2, and 
vary the e-ph coupling strength E^, as well as the band-filling. As for PBC, the largest system 
we can study for all fillings is the 14-site lattice, so we shall concentrate on this one at the 
beginning. 

In Figs. 12, we show the spectra for electron numbers = 2-7, phonon frequency 12 = 0.2, 
and sUghtly different coupUng strengths around = 0.6-0.75, i. e., values which are very 
close to the critical one at half filHng, E'^^iD. = 0.2) = 0.52. The common feature of all spectra 
is a dominant peak, which is slightly shifted away to higher co from the frequency (indicated 
by the dashed vertical lines) at which the non-interacting system shows its main response (see 
section III). At smaller frequency, there is also considerable absorption, which partly fills in 
the quasi-gap below the main peak. At low density (n^ = 2, 3), there is a rapid fall-off towards 
higher frequency, whereas for densities close to 1/2, we recover again the slowly falUng high- 
frequency tail, driven by the polaron-polaron interactions. At half-filUng, and for = 6, the 
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Figure 10: Schematic diagram of the dominant excitations induced by an external electromagnetic field 
in an (un)doped CDW. The doping of the CDW leads to the formation of mobile domain walls. 

CDW-gap, A^Dw beyond which a slowly falling high-frequency tail appears (if the coupling 
is not too large). As we dope away from half-filling, the dc conductivity is reestablished, and 
the gap disappears. Spectral weight is being transferred from high (co > A^^^) to low energy, 
but the slowly falling high-frequency tail remains. As we approach the low-density limit, the 
effective interactions become less and less important. This is reflected by the disappearence 
of the high-frequency tail, and the ratio of the Drude weight to the total spectral weight, which 
comes close to 1. For small coupUng (such that there is no CDW at n = 1/2), the doping depen- 
dence of ai(o) is not very pronounced, since the polaron-polaron interactions are negUgible. 
We shall therefore concentrate on the interesting case of sufficiently strong coupling. 
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Figure 1 1 : Doping dependence of the Drude weight 
D for various values of the coupling Eb obtained for 
a lattice with = 14 using (A)PBC at fixed phonon 
frequency Q, = 0.2 



0.0 



0.00 0.25 0.50 0.75 1.00 



n 



28 



Let us now turn to the calculations of the optical conductivity. As can be seen from Eq. (6), 
the evaluation of this quantity involves the whole set of excited states of the full Hamiltonian. 
Consequently, having a good variational approximation of the true ground state (obtained by 
optimization of the energy within a subspace of the full Hilbert space) does not guarantee 
that the optical conductivity calculated within the same subspace is as accurate. It can not be 
excluded a priori that excited states |/) outside of this subspace have a considerable current 
matrix element {i\ j|0) with the ground state within this subspace. Obviously, such contribu- 
tions are lost by our calculation, and it seems necessary to estimate their importance. Within 
our approach, the leading corrections can be obtained by extending the local phonon basis to 
additionally include the first excited states |1);, |1); of the (un)displaced oscillator. This will 
be done in the next section. 

In the accompanying paper [25], the doping dependence of the optical conductivity is 
investigated in the presence of PBC. There, we stress the strong influence of the polaron- 
polaron interactions on (7(ct)), which were neglected in previous calculations of this quantity. 
The generic physical picture that emerges is the following: For coupling strength large enough 
to result in an ordered CDW-state at half-filling, the conductivity at zero frequency, /. e., the 
Drude weight, is proportional to the number of electrons (as for a usual metal) at low electron 
(hole) density n <C 1/2, while at electron densities close to n = 1/2 (|n - 1/2| <C 1/2), it is 
proportional to the number of domain walls, = 2|ng -NI2\ (see Fig. 11). 

In the context of a CDW-state, a domain wall can be described as follows: Let us charac- 
terize the ordered CDW-state by a periodic function n{Ri) = 1/2 + cos(;r/?/a)/2 which mimics 
the alternating electron occupation on neighbouring sites. In this picture, a domain wall is a 
defect in the periodic structure, at which the phase (j) - nRila suffers a shift of K. To the right 
of this defect, all the elections are shifted by one lattice site (see Fig. 10). If an additional 
election (or hole) is intioduced into the CDW, it is easy to convince oneself, that this generates 
precisely two domain walls (at least in the presence of PBC), and they should behave hole- 
Hke, /. e., carry positive charge, for n < 1/2, and election-like, /. e., carry negative charge, for 
n > 1/2. Note, that the origin of this behaviour lies entirely in the effective electron-electron 
interactions mediated by the phonons. 

The doping dependence of (7((0) in this coupUng regime also reflects these interactions: At 
n = 1/2, the dc conductivity vanishes, and at finite frequencies o{co) shows the characteristic 
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Table II: The dimension of the total Hilbert space /?totW n, Q) as a function of the lattice size A^, and the allowed 
number of independent phonon states in each mode Q. 



Q 


N = 6 


N=S 


iV= 10 


iV= 12 


iV= 14 


N=16 


2 


1.280- 10^ 


1.792- 10^ 


2.580-10^ 


3.785-10^ 


5.623-10^ 


8.434-10^ 


4 


8.19210'^ 
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states of each mode. A detailed analysis of the results within these three subspaces will be 
published in a separate paper [25]. There, we use several different criteria to test the quality 
of each Ansatz: (i) The total energy, (ii) the ability to reproduce the exact parameters of the 
effective Hamiltonian in the strong coupling limit, and (iii) a comparison with available results 
from quantum Monte Carlo simulations. The result of this analysis is that in the interesting 
parameter range from intermediate to strong e-ph coupUng, there is a distinct subspace, which 
gives by far the best results. It is spanned by the direct product of the displaced (coherent 
state) 

\0)i=e-^"'^e^-i\0)i, (28) 

and the undisplaced oscillator ground state |0) ; at each lattice site. We refer to this variational 
subspace as Ansatz 11. 

In the calculations presented in this paper, we always used Ansatz 11. In the strong coupling 
limit, t/E^ 0, this choice guarantees the correct ground state energy, since the corresponding 
subspace contains the true degenerate ground state manifold which is obtained for t = 0. 
One has to be a Uttle careful with this choice, since the two mentioned basis states are non- 
orthogonal 

/(0|0); = e-^"'^ := V. (29) 

For the Lanczos-method, it is simpler to deal with orthonormal basis states, hence we choose 
the combination 

l<;>i )/ = - y\o)i ) ; 102)/ = 10)/ . (30) 

^ 1 _ y2 \ / 
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Unfortunately, the Lanczos method is not applicable to an e-ph system without further ap- 
proximation. The reason is simply the fact that the phonon Hilbert space is infinite-dimensional 
already for a single phonon mode, whereas the Lanczos algorithm can only be used for finite- 
dimensional Hilbert spaces, since convergence to the true ground state is only guaranteed if 
the vectors |4>^) = Y,m'^nm\^m) ^® expressed in a complete set of basis states of the full 
(in each symmetry sector) Hilbert space. This requires computer storage of the coefficients 

of at least two vectors at one time. For an infinite-dimensional Hilbert space, clearly, one 
would require an infinite amount of computer memory in order to obtain exact results. 

Therefore, the only possible way to apply the Lanczos method to an e-ph Hamiltonian 
is to use a variational approximation for the phonon Hilbert space, which makes it finite- 
dimensional. Since variational approximations are most often motivated by physical intuition 
or/and computational convenience, obviously, there are many possible choices for a variational 
ansatz. For our Hamiltonian, we considered the three 'most natural' variational approximations 
to the phonon Hilbert space, which are within the limits set by current computer power and 
the Lanczos algorithm. 

The basic computational limitation we encounter is the maximum number of independent 
phonon states that can be allowed for each mode. If we denote this number by Q, then on an N- 
site lattice, the total number of Unearly independent phonon states is given by Rp^iN, Q) = . 
On the other hand, the size of the electron Hilbert space (neglecting spin) on a A^-site lattice 
with n electrons is given by R^iN, n) = (^) , so that the total Hilbert space has the dimension 
R,ot(N, n,Q) - Q = . As an illustration, in table I, we show this number for different 

n\(N-ny. 

values of A^, Q, and fixed n - N/2, where the electron Hilbert space is maximal for a given N. 
Bearing in mind the size of core memory of present supercomputers (~ 1GByte, corresponding 
to a maximum size /?totW 2) ~ 4.5 • 10^), we can extract from these numbers that the 
maximum tractable lattice sizes are = 15 for g = 2, iV = 9 for g = 4, = 7 for Q = 6, 
and A^ = 6 for Q = 8. We have performed calculations mainly for Q = 2. The quaUty of this 
approximation, /. e., the influence of additional phonon degrees of freedom, was then tested 
by a comparison with results for Q- 4-. As we shall see later, the qualitative features obtained 
from the Q-2 calculation are not altered. The lattice sizes for Q > 4 that are tractable are too 
small to extract meaningful results. Therefore we did not consider them. 

Restricting ourselves to the case Q = 2, there are three plausible choices for the two phonon 
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tridiagonal matrix in the basis {|^„)}. In more detail, one obtains 



H\%) = ao\%)+l5o\^,), (24) 



The nth iteration is 

= Pn-l\^n-l) + + i^nl^n+l). (25) 



In practice, these operations are performed on a computer The main difficulty is to find an 
effective way for obtaining from |4>„). 

The extreme eigenvalues and -vectors of the such generated M x M tridiagonal matrix 
converge extremely fast to the real ones of H. Starting from a completely random |4>q) (which 
is the safest choice), one typically needs M = 60 - 100 iterations to obtain an accuracy to 12 
decimal places in the lowest eigenvalue Eq, and an error A = {{^q\H-Eq\^q))^''^ = 10"^- 10"^ 
for the ground state |^o)- The number of iterations needed for convergence of the ground 
state energy depends strongly on the splitting - Eq to the first excited state, and increases 
significantly as the two levels get close. 

For the calculation of the optical conductivity at zero temperature, one usually employs 
the continued fraction expansion of the Greens function [27] 

G(Z) = (Woly(Z-//)-Vl%) = ^"^'^^^^"^-2 ' (26) 

Z-Oo 



Z-a, - 



1 Z-a,--- 

where the coefficients (6c„, j8„) of the continued fraction are the usual Lanczos coefficients 
obtained from a second Lanczos run using the starting vector \^q) = JI^q). The continued 
fraction is truncated after M iterations (typically M = 300), when a good convergence of the 
spectrum is obtained. The real part of the optical conductivity is then given by 

ol^\(0) = — ^-Im G(<o + Er. + i5), (27) 
nco 

where we introduced a small imaginary part id, which gives a finite Lorentzian width to the 
delta peaks appearing in the original definition (6). 
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y. Numerical treatment of the t-J Holstein Hamiltonian 



A numerical method suitable for dealing with strong correlations would in principle be the 
quantum Monte Carlo (QMC) method. This technique was used by Hirsch and Fradkin [26] to 
investigate the CDW transition at half filling. They presented strong evidence that for arbritrary 
finite values of the phonon frequency Q., the spinless Holstein model exhibits a transition from 
a metallic to a CDW state at a finite critical value of the e-ph coupling constant ^g. They 
derived the t-V model as the effective second order Hamiltonian in the strong coupUng limit, 
and showed that the ratio of the effective parameters lim^^^^ tlV = 0. We recall, that most 
standard analytical approaches only take into account the hopping renormaUzation from the 
first order correction, which defines a small polaron band with a bandwidth ~ t = te~^^'^. As 
Hirsch and Fradkin noted, the fact that the t-V model exhibits a CDW transition at tlV - 1/2, 
gives strong evidence for an analogous CDW ground state in the strong coupling limit of the 
Holstein model. 

The drawback of the QMC method is that it is not straightforward to calculate dynamical 
quantities such as <7(<o), because the simulation gives direct results for dynamical quantities 
only as a function of imaginary frequencies. Therefore one needs to analytically continue to the 
real axis, which is not always a well-defined procedure. In contrast, the Lanczos method has 
been successfully appUed to calculate o{a)) at zero temperature for purely electronic Hamilto- 
nians on finite lattices [27]. This algorithm allows for the exact determination of the low-lying 
eigenvalues and -vectors of the Hamiltonian matrix. From the calculated ground state, it is 
then possible to extract dynamical quantities such as, e. g., the spectral function A(k, ft)), or 
the optical conductivity (7 (ft)), in principle with arbitrary accuracy. The big advantage of the 
method is that the results are exact, the disadvantage is the limitation (imposed by the com- 
puter memory available) to rather small lattices, so that often an accurate finite-size scaUng is 
impossible. 

The algorithm works as follows: One generates an arbitrary initial vector |4>q), with the 
only restriction that |^>q) should have a finite overlap with the ground state in the symmetry sec- 
tor of the Hilbert space one is interested in. Starting from |4>q), one generates an orthonormal 
set of states {!<!>„)} by successively operating with H on |<I>„), n = 0, • • • , M, and orthonormal- 
izing to all the previous vectors, such that the Hamiltonian is represented by a real symmetric 
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well as in the strong coupling case [24] using perturbative methods. However, in none of the 
approaches used were the polaron-polaron interactions taken into account, i. e., the results are 
only trustworthy (especially in the strong coupling limit) in the single polaron limit. 

There are a number of reasons, why we consider these approaches do not to apply for our 
problem: (i) We are interested in a situation with high carrier concentrations, e., polaron- 
polaron interactions are important, particularly in the vicinity close to density n - 1/2. (ii) The 
estimated value of the e-ph coupUng corresponds neither to weak nor strong coupUng, so there 
is no guarantee that perturbation series rapidly converge, (iii) We are also interested in the 
effect of the spin degrees of freedom, even though we do not expect them to be very important 
at reaUstic values of the spin-exchange constant, (iv) We would like to include the effect of 
vacancy disorder by choosing open boundary conditions as already explained in section III. 

An analytical treatment in the interesting parameter regime including the additional ef- 
fects of spin-exchange coupling and disorder is hopelessly difficult. Consequently, we chose 
to study the model using a numerical technique, the Lanczos method. In an accompanying 
paper [25], we present a detailed discussion of our numerical results for the simple Holstein 
model and contrast them with the analytic results mentioned above. 
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We have chosen the coupling to the hole, because in the t-J model, the unoccupied site cor- 
responds to the Cu™ singlet for which the coupling should be much stronger than for the 
Cu^ configuration which corresponds to an occupied site. The contribution to the coupling 
constant X from this mechanism can be obtained from the slope of the curve in Fig. 9 as 
X :- <iAj,ijjg|gt(Q)/(ig|g^Q = -4.2eV/A, where we defined Q to be positive for displacements of 
the 0(4) towards the Cu(l) atoms. However, this value does not contain the purely electro- 
static interaction energy, which should have opposite sign. We can estimate it to be of the 
order of 1 - 2eV/A. A rough value for the total X can thus be given as A = 2 - 3eV/A. 
The phonon part of the Hamiltonian is the usual one (neglecting dispersion) 

(21) 

with the reduced oxygen mass M = MqJ2 ~ 8amu, and the spring constant K. Li terms of the 
phonon creation and aniuhilation operators a J, a; 

with the phonon frequency Q. = \JkIM, the total Hamiltonian becomes 

H = Hj_j +na a]ai + A ^(1 - n;,)(a; + a]), (23) 

/ Is 

and the new coupling constant is A = (JiX^ I2MQ)^'^ . The phonon frequency for the mode under 
consideration is known rather accurately from Raman scattering [19, 20] as/iQ ~ 60-65meV, 
depending on the material. This allows for an estimate of the spring constant as ^ = Sl^M » 
7-8eV/A^. The polaron binding energy := X^/{2K) = A^/{hQ), which is the most convenient 
way to measure the effective interaction strength, is therefore of the order of = 0.25-0.64eV. 
This is a value comparable to the effective bandwidth 4t ~ 1.6eV, and therefore puts us in a 
regime of intermediate coupling, which is the most difficult to treat analytically. Note, that our 
choice of coupling the phonon displacement to the local hole rather than electron density does 
not affect any properties of the e-ph Hamiltonian (23). Apart from a constant energy shift, the 
Hamiltonian is electron-hole symmetric (as long as the kinetic energy part is). 

The spinless version of the Hamiltonian (23) is the usual Holstein or molecular-crystal 
model [21]. The optical conductivity of this model has been studied in the weak [22, 23] as 
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Figure 8: Schematic diagram of (a) the vibrational mode which involves displacement of the apical 
oxygen in the c-direction towards the neighbouring Cu(l)-site, and (b) the coupling of this mode to the 
Cu™ singlet. 

From the plot of I^^^^^J^Q) := ^singietCG) - ^singiet(2 = 0) in Fig. 9, we note that for values 
Q < 0. 1 A, the curve is pretty linear, and only for larger Q, the curvature becomes important. In 
view of this, it seems appropriate to model the e-ph interaction by an additional Holstein-like 
term in the Hamiltonian, which couples the local displacement Qi linearly to the local hole 

density 1 - 

^e-ph = AE2/(l-'^/.)- (20) 
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Figure 9: The change Agmgiet in total energy of the 
Cu™ singlet state as a function of the displacement Q 
of the apical oxygen as shown in Fig. 8. The gain in 
kinetic energy of the singlet leads to an estimate of the 
e-ph coupling constant of A = -4.2eV/A. 
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IV. Polaronic Electron-Phonon Coupling 



It is well known from various studies of small polaron physics, that systems with a strongly 
localized electron-phonon (e-ph) interaction usually show a infrared absorption band. There- 
fore it is worthwhile to look more closely at possible phonon modes in the 1-2-3 materials, 
which might give rise to such a coupling. The standard theoretical approaches to these sys- 
tems, employ some kind of mean-field decoupling to handle the polaron-polaron interactions. 
This can give satisfactory results in the limit of small carrier density, or coupling. However, in 
the case of interest, the carrier density n ~ 0.5/Cu is not low, and even close to a commensurate 
value at which it can be favourable to form some kind of a lattice superstructure, i. e., a charge 
density wave (CDW). Furthermore, as we shall show below, the estimated strength of the e-ph 
coupling in the CuOg -chains puts us in a regime of intermediate coupling, which is difficult to 
handle with standard perturbative methods. Hence, we need a technique, which is able to take 
into account the effect of such correlations in a satisfactory way. 

We start from the ID t- J -model as the effective electronic model having been motivated 
above. Bearing in mind that an empty site in the t-J-mode\ corresponds to a Cu"^ singlet 
formed from Cu 3d^2_y2 and O 2py^ holes, we can identify an effective way of coupUng the 
apical oxygen 0(4) vibrational mode with displacement in the c-direction to the holes in the 
singlet (see Fig. 8): The energy of the singlet is mainly governed by the two hopping matrix 
elements tp^ between Cu 3d^2_y2 and O 2py ^ orbitals, and tpp between neighbouring O 2p^^ 
orbitals. As the two apical O atoms above and below a Cu are moved towards it, the values of 
tpj, tpp increase, and the total energy of the singlet is lowered due to the gain in kinetic energy. 
Taking for the singlet the approximate wave function 

|Cu™) = Uq cl^cl^ + ai [cl^cl^ - cl^cl^^ + ^icl^cl^ \), (18) 

which is localized on a single Cu and its four surrounding oxygens, we can calculate the 
total energy of the singlet ^singlet ^ function of the displacement Q (measured from the 
experimental Cu(l)-0(4) distance) using the empirical scaling laws for the hopping matrix 
elements [18] 

^pd(^pd> _ I ^ \ . hp^^pp^ _ ( ^'pp \ 

^pd^^pd) ypdj ^PP^^PP^ VPP J 

where r^^, r^^ are the relevant interatomic distances between atoms with dor p orbitals. 
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wit < 2, and > 8, we recover in the high frequency tail, the structure of the isolated quasi 
Drude peaks which correspond to very short chain segments. 

A comparison of e. the = 8 calculation to the Pr spectrum, shows that now only a 
value ofW~ 1.2eV is needed to obtain a fit as good as the one for = 0, W/t = 4. However, 
going to an even larger n^- 16, the correct peak position requires a value of W/t = 1, in which 
case the high-frequency tail is not very well fitted anymore, especially regarding the strong 
quasi Drude peaks in the calculated spectrum. 

In summary, a model including both potential and vacancy disorder exhibits an optical 
conductivity spectrum which is extremely close to the measured one. However, the potential 
disorder strength needed to fit the experiment, W~1.2- 1.6eV, is rather large, and the source 
of such a strong random potential on each site is not obvious. A possible type of defect 
producing such a strong potential could be an oxygen interstitial between two Cu(l) belonging 
to neighbouring chains (see Fig. 1). The corresponding electrostatic potential energy = 
e /(e r) can be estimated to be ~ 1.2eV (using a rather small value for the dielectric constant 
e = 5) at the regular oxygen sites nearest neighbour to the interstitial. Nevertheless, it seems a 
rather high density of such defects would be necessary to produce the required random potential 
on each site. 

In the next chapter, we shall investigate a third mechanism which could give rise to the 
observed optical response, namely polaronic electron-phonon coupling. 
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Figure 7: Optical conductivity from the model with potential and vacancy disorder. The results were 
obtained for a chain of length A'^ = 100 at n = 1/2, using different values for the disorder strength W, and 
vacancy number rid- 



6, 7, show the spectra for defect numbers - 4, 8, 12, 16, and different strengths of potential 
disorder W. If we compare different values of for fixed W, we observe that the peak shifts 
to larger energies as is increased (as expected). At the same time, for smaller values of 
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Figure 6: Optical conductivity from the model with potential and vacancy disorder. The results were 
obtained for a chain of length A'^ = 100 at n = 1/2, using different values for the disorder strength W, and 
vacancy number rid- 



fit the experiment. We shall comment on this point at the end of this section. 

Let us now proceed by including vacancy disorder, in this case restricting ourselves to 
values of ^It = 0.1, V = 0, since the effect of a finite V was found to be not so crucial. Figs. 
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Figure 4: Optical conductivity from the model for pure potential disorder. The results were obtained for 
a chain of length N =100 atn = 1/2, and different values for the disorder strength W. 
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Figure 5: A comparison of the experimental spec- 
trum obtained in Ref. [1] with the curve calculated 
using the model with pure potential disorder, and dis- 
order strength W/t = 4. 
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CO ~ t/N of such segments shifts to higher energy. At the same time, the spectra become more 
and more noisy, dominated by strong isolated peaks corresponding to certain values of A^. 

Setting t = 0.4eV, we can try a comparison with the experiment. The present model 
predicts a shift of the observed mid-infrared peak to higher frequencies, when the oxygen con- 
centration is reduced. Since the exact oxygen content (and therefore n^) is difficult to estimate 
experimentally, we cannot say with certainty whether the absolute value of the calculated peak 
position is at the correct frequency for a particular sample. However, assuming an upper bound 
for the oxygen deficiency 5 < 0.1 in REBa2Cu307, the calculated value appears at ~ O.leV, 
half of the experimental value. In addition, there is also no sign of a smooth high-frequency 
tail as observed experimentally. Furthermore, the data in reference [3] taken from samples 
YBa2Cu307_5 with different 5 do not show any strong influence of 5 on the peak position, in 
contradiction to the calculations. In conclusion, it seems that the oxygen vacancy disorder (at 
least if modeled as we did) is also not enough to explain the observed mid-infrared peak. 

The second kind of disorder that could be present is simple potential disorder, which might 
originate, e. g., in randomly positioned oxygen interstitials (see Fig. 1). We model this by 
random on-site energies on each site, uniformly distributed in some interval [-W/2, W/2], 
the size of which is a measure of the disorder strength. In addition, we can also allow for 
variations in the hopping matrix elements as before. This type of model is often used in the 
context of Anderson locaUzation, and it is known from various studies [17], that a{co) shows 
a peak at low energy, and a high-frequency Drude-hke tail. Nevertheless, we have calculated 
g{(0) also for such disorder, since we were interested in the influence of the random f;. 

We first study the effect of pure potential disorder, setting - 0. Fig. 4 shows the spectra 
calculated for N= 100, ng = 50, and various values of W/t. As expected, we find a peak at low 
energy, and a subsequent slowly falUng high-frequency tail. The position of the peak rapidly 
shifts to higher frequencies as the disorder strength W is increased. Setting t = 0.4eV, we find 
that a value ofW~ 1.6eV is needed, to obtain the experimental peak position of ~ 0.2eV. 
In Fig. 5, we superpose the spectrum from PrBa2Cu307_5 with the calculated one at W/t = 4. 
The agreement in both, the shape and position of the peak, as well as the high-frequency tail is 
astonishing. From a purely theoretical point of view, it therefore seems that this model is fully 
adequate to describe the optical conductivity of the CUO3 -chains. However, doubts concerning 
the vaUdity of this scenario arise because of the large disorder strength W ~ 1.6eV required to 
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Figure 3: Optical conductivity from the model for vacancy disorder. The results were obtained for a 
chain of length A'^ = 100 at n = 1/2 using the parameter t'/t = 0. 1, and different values for V, and n^. 



the vacancies, and therefore effectively cuts the chain into isolated pieces. 

As we increase the number of defects n^, we observe that the low-energy peak broadens, 
and shifts to higher frequencies. This can be explained from the fact that the average length 
of perfect chain segments decreases as grows, and consequently, the quasi Drude peak at 
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Figure 2: Optical conductivity from the model for vacancy disorder. The results were obtained for a 
chain of length A'^ = 100 at n = 1/2 using the parameter t'/t = 0. 1, and different values for V, and n^. 



The characteristic feature of all spectra is a broad peak at low energy, and a rather rapid 
fall-off towards higher frequencies. Comparing the results for different values of V with fixed 
n^, t'/t, we note that the influence of V is not very big, when varied between V/t = 0- 8. The 
relevant effect seems to result from the reduced t' which introduces strong back-scattering at 
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with 

= ^h iW+l,s - (16) 

With the help of Eq. (15) it becomes very easy to calculate oi(o) at T = 0, since creates only 
single particle-hole excitations above the ground state, and hence (7[*^(<o) can be written as 

2 I ~ |2 

cjT\(o)^^ J2 j^5[co-(e,-l)]. (17) 

i V Cj. 

s < ne,r > He 

In all of the calculations of this section, we shall use open boundary conditions. Consequently, 
the Drude weight must vanish (more precisely, it is shifted to a delta peak at a frequency of 
the order t/N), and the entire spectral weight appears at finite frequencies. However, for chain 
lengths > 100, one has essentially reached the thermodynamic limit for the shape of ai^{(0) 
at frequencies co <c t/N, in which we are interested here. 

A realistic model for a situation where the oxygen vacancies are the dominant source of 
disorder (as seems to be the case for the CuOg -chains) can be constructed by an appropriate 
choice of the coefficients e^, tf. For a single vacancy between Cu sites / and 1+ l,w& assume 
a strongly reduced hopping matrix element t'/t <C 1, and in addition, a screened Coulomb 
potential which we can model, e. by a repulsive e^, e^+j = V > possibly extending to the 
neighbouring sites 1+2, l-l, etc. with an exponential tail. According to the oxygen deficiency 
in REBa2Cu307_5, we can then assume a certain concentration 5 = nJN of such defects (n^ 
the total number of defects), and distribute them randomly. 

In Figs. 2, 3, we show the spectra for a chain of length N - 100 with t'/t -0.1, and different 
values of V, and n^. A length of N = 100 is sufficient, since we are not interested in the very 
low frequency behaviour. We compared the results for chain lengths of 100, 200, and 400, 
and found no significant difference at frequencies co > t/50. The precise value of t'/t does not 
affect the shape of the spectra as long as t'/t <C 1 (which is the physically relevant region). In 
these calculations, we used a potential which extends to the three sites neighbouring a vacancy 
falling off exponentially with a screening length equal to the lattice constant. 

We always used a band filling of n :- nJN - 0.5, in accordance with the experimentally 
estimated value mentioned in the introduction. Each curve represents the average over 1000 
samples of the disorder, and for clarity of presentation, we introduced a smoothening of the 
delta functions by averaging the points over small frequency intervals typically of size t/200. 
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Figure 1 : Schematic diagram of an oxygen vacancy (a) and an oxygen interstitial (b) in the CuOs-chains. 
/fjjjg is diagonal in the single-particle states 

r 

r r 

The coefficients are 

i0/)=cjio) i0,)=c;io). 

The many -particle eigenstates are simple Slater determinants 

i*(K}))=n^«jo)' 

{nr} 

where {n^} specifies single particle levels {n^ = number of electrons). If we label the single- 
particle levels in ascending order (e^ < e^^j), the ground state is given by 

1%) = n ^r\0)- (14) 

r<ne 

In order to calculate expectation values of j^, we express in terms of the 
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III. Disorder Effects 



Various experimental studies of the high-r^ oxides in general, and the 1-2-3 materials in par- 
ticular have shown that it is very difficult to control the oxidation, /. e., the oxygen content of 
a sample. The high mobility of the oxygen atoms (in some compounds even down to below 
room temperature) makes it practically impossible to obtain exactly stochiometric compounds. 
This is especially severe in the case of the REBa2Cu307_g family, where the 0(1) atoms which 
form the Cu-0 bonds along the Cu03-chains, have an extremely high mobility. Therefore, one 
expects (and finds) a large degree of disorder in the chains, in particular originating from O 
vacancies, and possibly interstitial s. 

A first extension of the f-J-model should therefore be the incorporation of disorder. Since 
we have seen in the last section that in ID for reaUstic values of Jit ~ 1/4 the spin exchange scat- 
tering hardly alters the optical conductivity as compared to non-interacting spinless fermions 
(corresponding to / = 0), we will neglect the /-term in the calculations of this section, i. e., 
we treat the electrons as spinless. This simpUfies the calculations dramatically, because the 
problem reduces to a single-particle one. 

We will consider two kinds of disorder: (i) disorder in the hopping matrix element t, and 
(ii) randomly disordered on-site energies Type (i) should be a realistic model for oxygen 
vacancies (see Fig. 1) taking into account that across a vacant oxygen site the effective overlap 
matrix element will be drastically reduced. Furthermore, a vacancy can also be expected to 
produce a strong electrostatic potential arising from the missing charge of the O"'^ ion. A 
possible source of potential disorder leading to random on-site energies e^, are, e. g., oxygen 
interstitials (see also Fig. 1) or other defects. 

The spinless Hamiltonian including the effect of disorder can be written as 



where the ti, are random variables with certain probability distributions which will be spec- 
ified later. This Hamiltonian contains only single-particle terms, and therefore its spectrum as 
well as the associated single-particle orbitals are easily obtained by numerical diagonalization 
of the corresponding tridiagonal matrix. For the calculation of al^^((o), we can still use Eq. (6) 
with the modified current operator 




(9) 




(10) 
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oico) is dominated (> 99% of the total weight) by a Drade peak, and only very little weight 
at finite a>. Of course, as we cross the phase separation line for very large values of / > 3t, 
the situation changes: (i) The total weight decreases with increasing /, and (ii) more and more 
spectral weight shifts from the Drude peak to finite (O, until we end up with an insulator at 
J - oo. But these values are far too large to explain the experimentally observed spectrum. 
We therefore conclude that the ID t- J -model alone cannot give a satisfactory explanation for 
the optical response of the Cu03-chains. So we have to include additional terms in the basic 
Hamiltonian (2). 
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where ^q, ^„ are the eigenstates of the Hamiltonian (2) at ^> = 4>o. For the total integrated 
spectral weight a^^^ under o^{a>), there exists an /-sum rule 



with T the kinetic energy in (2). With the help of this sum rule it is possible to give an 
alternative expression for the Drude weight D 



which is helpful to check the calculated results. 

This formalism has been used to calculate c7(<o) for the ID t-J-model [7] numerically on 
finite chains employing the widely used Lanczos algorithm. The (at first sight) surprising result 
was that the scattering introduced by the strong correlations is barely noticeable: The integrated 
spectral weight at finite co appears to be an extremely small fraction of the total weight o^Q^, 
which indicates that the optical response in the ID t-J-model is almost identical to the case of 
spinless fermions. The deeper reason for this is the separation of spin and charge, which is a 
well-known and generic feature of ID models with interacting electrons. These models do not 
exhibit Fermi-Uquid behaviour, but belong to a different universaUty class, so-called Luttinger 
liquids [14]. 

For the ID f-/-model, an exact solution exists only at the supersymmetric point / = 2t, at 
which the difference in the charge and spin velocity can be rigorously shown [15]. However, in 
the small J/t Umit, we approach the strong coupling limit of the ID Hubbard model, for which 
only quite recently Ogata and Shiba have shown [16] that the ground state Bethe-ansatz wave 
function factorizes into a charge and spin part. The charge part corresponds to the ground state 
of the corresponding system of spinless fermions, and the spin part to the ID AF Heisenberg 
ground state. In fact, the spin-charge decoupling can be intuitively understood from the fact 
that in ID, the charge motion alone is not capable of changing the spin configurations, since 
this requires the existence of closed loops in the lattice topology. 

This factorization, e., the charge-spin separation, has the effect, that in the small J/t limit, 
the optical conductivity consists basically of the one expected from non-interacting spinless 
fermions, since the electric field couples only to the charge degrees of freedom. Consequently, 




(7) 




(8) 
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In this effective description, / labels the Cu(l) sites of the chain, cj^ creates a hole in the Cu 
3<i^2_y2 orbital with spin s, its neighbouring O 2p-shells being filled (Cu°), and C;^ creates 
a singlet formed from Cu 3d^2_y2 and O 2/7^^ holes (Cu™). The projection operators "P, = 
l-ni^rij^ eliminate configurations with doubly occupied sites. Again, we argue that the effective 
parameters t and / should be very similar to the ones derived for the plane [10], and hence, 
we adopt those values, t ~ 0.4eV, J ~0. leV. 

We use the well-known Kubo formula for the calculation of the optical conductivity a{co). 
In particular, for a one-band tight-binding Hamiltonian like (2), the prescription to calculate 
g{co) at wave vector q = is the following [13]: If we use periodic boundary conditions 
(PBC), the real (absorbent) part (7i(ct)) of a t-J ring can be decomposed into a Drude and a 
regular term 

(7i((») = Ddia) + Oi\q)), (3) 
where the Drude weight D is given by 



D=kN 



(4) 

3>=3>o 



with Eq{^) being the ground state energy of the system in the presence of a time dependent 
magnetic flux 4> through the ring, and N is the number of lattice sites. The flux can be generated 
by a uniform vector potential along the ring (x-direction), A = Ae~'"'e^. (4> = NAel{hc)) which 
modifies the hopping by the usual Peierls phase factor t r\ tGxp{±i^/L). The lattice constant 
Uq has been set to 1. <I>q is the flux for which the total energy is minimal. Usually either 
4>Q = (corresponding to PBC) or 4>q = ;r [corresponding to antiperiodic boundary conditions, 
(APBC)], depending on the band filUng. For the ?-/-model, e. g., one obtains 4>q = if the 
number of electrons satisfies = Am + 2, but ^q- it for = Am (m a positive integer). In 
general, at 4> = 4>q, the expectation value ( y^) of the current operator at q = 

l,s 

in the ground state |^o(^o)) vanishes, which is required to make (4) appUcable. The regular 
term a[^^((u) can be written as 

^1 ^ ^ 2^ — -B 'B ~ " -^o)]' (6) 

■'V — ; ii™ — /in 



n. The Electronic Structure of CuOg -chains 



When trying to construct a model for the electronic structure of the CUO3 -chains, the first thing 
one should note is the fact that the chain Cu(l) site has the same four-fold oxygen coordination 
as the planar Cu(2,3) (neglecting the buckling of the 0(2,3) atoms out of the plane). Therefore, 
we can expect a similar Cu-0 bonding pattern in both cases. This is confirmed by band- 
structure calculations for YBa2Cu307_5 [8], which apart from the common strongly dispersive 
planar CuO pd<7-hand also find a chain-band crossing the Fermi energy Ep. This band is 
formed from an antibonding pdo combination of 0(l)2p^-Cu(l)3(i_2_,.2-0(4)2p^ orbitals, and 
shows a dispersion similar to the planar band, but only along the chain axis. 

Regarding the similarity between planes and chains, it seems reasonable to try a description 
of the electroiuc structure by the ID analog of the Emery model [9] 

ijs i 

where i,j label the Cu(l) and 0(1,4) sites, the operators c]^ create ?)d^2_y2-C\y and 2py^-0 holes, 
and s is the spin index. The include the usual on-site energies £j - (by convention), 
= 3.6 (which we assume to be the same for 0(1,4)), and the nearest neighbour (n.n) hopping 
integrals t^^ = 1.3, t^^ = 0.65 (all values in eV). The cited values correspond to the ones found 
for the plane [10], which, regarding the above mentioned similarities, and the similar CuO 
bond lengths, also seem reaUstic for the chains. Further justification for the use of these values 
comes from measurements of the optical gap in the compound CajCuOg [11], in which one 
finds the same CUO3 -chains with a Cu^ oxidation state, i. e., in the insulating charge-transfer 
regime. The observed value of the optical gap is ~ 2eV, which agrees well with the values 
found in the planar antiferromagnetic charge-transfer insulators like YBa2Cu306, and therefore 
suggests a similar value of - e^. 

In the Ught of this, it seems appropriate to apply the same procedure, as for the CuOj- 
planes [12], to derive an effective low-energy Hamiltonian from (1), which consists of a single, 
strongly correlated tight-binding band (projected on the subspace of singly occupied sites), and 
a nearest-neighbour antiferromagnetic exchange interaction, /. e., the ID t-J-model 

l,s l,s ^ ^ 
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rather strongly coupled to the electrons. Using the t- J -model as the effective electronic model, 
we will show that this coupling can be included in the form of a Holstein term, in which 
the local ionic displacement couples linearly to the local charge. Employing a variational 
approximation for the phonon Hilbert space, we calculate (7(g)) for finite systems using the 
Lanczos technique. The effect of the vacancy disorder can be incorporated by the choice of 
open boundary conditions. 

Finally, we will discuss the results of our different approaches in the light of the experi- 
mental data, and also suggest possible experimental tests of the different mechanisms. 
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tail falls off much slower than the simple Drade formula {p ~ II (O at high frequency) would 
suggest. At energies below 0.2eV, there seems to be a difference between the two compounds: 
while in YBa2Cu307_g, the conductivity joins onto a Drude-peak above T^} in PrBa2Cu307_g 
it drops sharply to zero, reflecting the insulating nature of the latter compound. Note, that in 
the Pr-compound, this behaviour of the conductivity has important implications for possible 
mechanisms explaining the observed rp-suppression[6]. The above described shape of Oc^((») 
clearly indicates that the electronic structure of the chains is non-trivial, and it stimulated our 
interest to find a possible explanation. 

In this article, we offer two likely mechanisms which could lead to the observed o^^^G)): (i) 
the presence of very strong potential disorder, and (ii) moderately strong polaronic electron- 
phonon (e-ph) interactions, possibly combined with the effect of vacancy disorder While 
the first explanation gives a better (in fact excellent) fit to the experimental data, the second 
one seems to be more reasonable from a physical point of view, since, as we shall show, the 
potential disorder strength needed to obtain agreement with the experimental data is hard to 
justify. 

The plan of this paper is as follows: First we shall identify the one-dimensional (ID) t-J- 
model to be a good starting point for the description of the electronic structure of the CuOg- 
chains. However, as is well known by now from various studies, the optical conductivity of 
the ID ?-/-model is more or less identical to that of spinless fermions (as long as / is below the 
threshold for phase separation) because of charge-spin separation. There is barely any spectral 
weight at finite frequencies [7]. Thus, one is forced to include other interaction terms in the 
Hamiltonian, which could possibly reproduce the experimentally observed strong scattering. 

We continue with a simple extension of the model by including the disorder from oxygen 
vacancies on the bridging 0(1) site. This is expected to be the main source of disorder in the 1- 
2-3 materials. It is mainly due to the large mobility of the oxygen atoms at higher temperatures. 
Such vacancies can be modeled by interrupted ID t-J chains, but we will show that such a 
scenario cannot reproduce the observed o^^{g)). However, as mentioned above, the inclusion 
of strong potential disorder turns out to yield an excellent fit to the experiment. 

As an alternative explanation, we then consider moderately strong e-ph interactions. We 
suggest that the apical oxygen mode with a displacement in the c-direction is likely to be 

^This point is still somewhat controversial. More recent reflectivity measurements of YBa2Cu307-5 down to lower 
frequencies [5], indicated that even above Tc there is no upturn in the conductivity as the dc limit is approached. 
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I. Introduction 



The 1-2-3 family of high-r^ oxides with its most prominent member YBa2Cu307_5 is the 
only example of the new superconductors which in addition to the common CUO2 planes also 
contains one-dimensional CUO3 -chains contributing to the electronic density of states at the 
Fermi level. ^ At the same time, this class of compounds is the best characterized among 
the high-r^ materials, and furthermore, many of the technological appUcations are based on 
YBajCugOy.^. However, in the interpretation of many experimental data, it is not obvious 
how the features of the chains can be separated from those of the planes. The latter are known 
to drive the superconductivity in all high-T^ compounds, and are therefore of primary interest. 
Only recently, with the availabiUty of single-domain crystals, a clearer picture of the interplay 
between chains and planes has emerged. In any case, for a complete understanding of the 
properties of the 1-2-3 compounds, it is absolutely necessary to distinguish, and where possible 
isolate the features of the CUO3 -chains. 

One of the physical quantities which contains valuable information about the electronic 
structure of a system is the optical conductivity a(co). Its frequency dependence is determined 
by both, the character of the ground state, as well as (in principle) all the excited states. 
Experimentally, (7(o)) can be determined quite easily, e. g., from reflectivity measurements 
by applying a Kramers-Kronig transformation. 

The synthesis of detwinned YBa2Cu307_g, and PrBa2Cu307_5 crystals made it possible 
to extract the anisotropy of (7(<o) in the crystallographic a- (perpendicular), and b-direction 
(parallel to the chains), by polarizing the electric field E in one or the other direction. This 
revealed an interesting and important property of these compounds, namely that a^(G)) (with 
E II b) shows a strongly enhanced spectral weight compared to cT^/co) especially in the infrared 
region ((o < 1.5eV). A procedure which suggests itself is then to subtract (y^i(o) from o),(<o), 
and to associate the corresponding oscillator strength with the CUO3 -chains, the main source 
of anisotropy in these materials. 

In aU pubUshed data (Ref.[l] for PrBa2Cu307_5, and Ref.[2, 3, 4] for YBa2Cu307_5), 

these difference spectra o^^ico) := (T^(cd) - C7^(©) show a characteristic broad mid-infrared peak 

centered around 0.2eV, and a slowly falling high-frequency tail extending beyond leV. This 

^In fact, this statement is not quite correct. There exists a slight variation of these materials, the so-called 1-2-4 
compounds, with its prototype YBa2Cu408. They contain one-dimensional CuOs double-chains, which also provide 
states at the Fermi level. 
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Abstract. The optical conductivity of the CuOs-chains, a subsystem of the 1-2-3 materials, is dominated 
by a broad peak in the mid-infrared (o) = 0.2eV), and a slowly falling high-frequency tail. The ID t- 
/-model is proposed as the relevant low-energy Hamiltonian describing the intrinsic electronic structure 
of the CuOs-chains. However, due to charge-spin decoupling, this model alone cannot reproduce the 
observed (7(eo). We consider two additional scattering mechanisms: (i) Disregarding the not so crucial spin 
degrees of freedom, the inclusion of strong potential disorder yields excellent agreement with experiment, 
but suffers from the unreasonable value of the disorder strength necessary for the fit. (ii) Moderately 
strong polaronic electron-phonon coupling to the mode involving Cu(l)-0(4) stretching, can be modeled 
within a ID Holstein Hamiltonian of spinless fermions. Using a variational approximation for the phonon 
Hilbert space, we diagonalize the Hamiltonian exactly on finite lattices. As a result of the experimental 
hole density « 1/2, the chains can exhibit strong charge-density-wave (CDW) correlations, driven by 
phonon-mediated polaron-polaron interactions. In the vicinity of half filling, charge motion is identified 
as arising from moving domain walls, i. e., defects in the CDW. Incorporating the effect of vacancy 
disorder by choosing open boundary conditions, good agreement with the experimental spectra is found. 
In particular, a high-frequency tail arises as a consequence of the polaron-polaron interactions. 



PACS Numbers: 71.38, 71.45.L, 74.25, 74.72, 78.30, 78.30.L 



New Address: Science and Technology Center for Superconductivity, Materials Science Division, Argonne National 
Laboratory, 9700 South Cass Avenue, Argorme, Illinois 60439 



